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Abstract 

Recently, the LHCb collaboration reported the first evidence for the decay B® — > fi + fi~. A 
branching ratio of B(B® — > u + u~) = (3. 2^2) x 1CP 9 is given. Using the newest data, and together 
with the most precise predictions of the Standard Model contributions to the decay, we derive 
the constraints on the combinations of the R-parity violating parameters. Our results are several 
orders of magnitudes stronger than the constraints in the previous literature. We also update the 
constraints on the relevant parameters using the upper limit of B(B^ — > fi + fi~). 

PACS numbers: 12.60.Jv, 13.20.He, 14.80.Ly 
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I. INTRODUCTION 



The helicity suppressed rare decay process B® — > fi + fi~ is induced by Z boson mediated 
penguin diagram and box diagram in the Standard Model (SM). This double suppression 
mechanism makes the SM prediction for the process very small 1| 

B(B° S = (3.23 ± 0.27) x 1(T 9 . (1) 

The fact that there are only leptons in the final states makes it a golden channel for the 
discovery and/or constraining the new physics model parameter space, since the new physics 
contributions can be larger than the SM effects and there is the least hadronic uncertainty. 

The minimal supersymmetric standard model with R-parity violation (MSSM-RPV) is 
an extension of the minimal supersymmetric standard model (MSSM) by abandoning the 
discrete symmetry, the R-parity, which is defined by R p = (^—i^ B + L + 2S ; where B is the 
baryon number, L is the lepton number, and S is the spin of the particle. The most gen- 
eral R-parity violating term can be included in the MSSM by introducing the following 
superpotential 

W K = HHvLi + \\ l3k UL 3 El + KjkhQjDt + ^\" jk U?DjDl, (2) 

where additional factor of 1/2 is added because of the fact that the first two indices of the 
couplings Xijk and A" jfc are antisymmetric. It is easy to see from this superpotential that the 
B° s — y decay can be induced at the tree level from the lepton number violating terms 

A and A'. 

Study of the MSSM-RPV have been performed in many rare decay processes. The bounds 
on the relevant parameters in the MSSM-RPV obtained from the decay B® — > fi + fi~ were 



derived in Ref. 



4j and revised in the literature However, there were only upper 



bounds from experiments at that time. 



Recently, the 



of B° a -+ n+ti- E] 



jHCb collaboration reported the first measurement of the branching ratio 



B(B° S fa-) = (3.2l}j) x 1(T 9 . (3) 

This just lies on the central regions of the SM prediction in Eq. ([T]), which will put severe 
constraints on every new physics models. In this brief report, we will use the newest data 
in Eq. ([3]) to constrain the relevant parameters in the framework of the minimal super- 
symmetric standard model with R-parity violation. Since the experimental measurement 



is quite close to the SM prediction, we have to include the contributions of the standard 
model together with the new physics contribution. Using the first time measurement of the 
branching ratio of the decay process, we give the most stringent constraints on the relevant 
parameters in the MSSM-RPV. We also update the constraint on the sneutrino exchange 
term and the squark exchange term from the newest experimental upper limits [7j: 

B(B° d -> < 9.4 x 1(T 10 . (4) 

This brief report is organized as follows, in Sec. HI] we present the analytical expressions; 
and then we use these equations to give the numerical results and discussions in Sec. IIII1 
We close this paper with a conclusion in Sec. IIVI 



II. ANALYTICAL EXPRESSIONS 



In this section, we present the formalism for the calculation of the branching ratio of the 
B° s — > in both the SM and the MSSM-RPV. The same formalism can also be applied 

to the process B d — > fi + fi~ with changing the corresponding Cabibbo-Kobayashi-Maskawa 
matrix elements. In the SM, the effective Hamiltonian governing B® — > can be written 

as js ] 

WgJ = -- 1 - a V t tV tsVY Y (x t )(h,P L s)mYPLfi) + h.c, (5) 
V2 27T sin <d w 

with Pl = (1 — 75)/2. The function Y is the famous Inami-Lim function calculated from 
the electroweak penguin and box diagrams 

„ , . x ( 4 - x 3x . \ 

While r/Y = 1.026 ± 0.006 includes higher order QCD corrections [lo| . 
The branching fraction of B Q S — > can be given by 



where t(B®) is the lifetime of the B° s meson. 

In the MSSM-RPV, the relevant effective Hamiltonian can be obtained by matching the 
amplitudes in full theory as shown in Fig. [T]onto the effective four fermion operators. For 
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FIG. 1: Feynman diagrams contributing to B® — > in the minimal super symmetric standard 

model with R-parity violation. 

the left diagram of Fig. dj the effective Hamiltonian from the exchange of sneutrino can be 
written as 

H^ff = ~(fiP R »)(bP L s), (8) 

where Pr = (1 + 7s)/2 and 

A = J2^. (9) 

For this new kind of contribution, the interference of the (flPnfi) density with the standard 
model (Ji^Plh) current leads to zero, so we can directly separate the SM and the MSSM- 
RPV contributions to the branching ratio as 

B(B° S -> // V) = & SM + Ba PV , (10) 
where we will use Eq. (CQ) as input value of B SM in our numerical calculations, while 



For the up-squark contribution shown in the right diagram of Fig. [H the effective Hamil- 



tonian is 

L ef/ 



^57 = +^(h^p R s)(firPL^ (12) 



where 

B = J-^™. (13) 



For this scalar quark contribution has the same structure of current as the standard model 
case, and the SM prediction lies in the central values of the experimental data, we can 
reasonably assume that the interference term of the SM and the MSSM-RPV is greatly 



TABLE I: Input parameters for B® and B® mesons used in numerical calculations. Uncertainties 
of these parameters are not considered for the reason that the uncertainties induced by these 
parameters on the relevant constraints on the couplings of the MSSM-RPV are far beyond the 
scope of the experimental data. 





T£(PS) 


/b (MeV) 


m B (GeV) 


B° 

^ s 


1.466 


227 


5.36677 


B° 


1.519 


190 


5.27958 



larger than the pure MSSM-RPV contributions, so we will approximately write the total 
contributions as 

B(B° /iV) = B SM + B%*, (14) 

where we also use the numerical value for B SM from Eq. (Q]); while the interference term 
can be written as: 



r(B° s )G F a 



4m?, 



16y 2^ sin kd w 



The B° s meson decay constant /#o shown in the above equations arises from the calculation 
of the hadronic matrix element, which is denned as 

^\h^s\B° s )=tf B ^ Br (16) 

The hadronic matrix elements for the pseudo-scalar density can be derived from the equation 
of motion under the assumption that m& ~ m^o. 



III. NUMERICAL RESULTS 



In this section we present our numerical results. Following Ref. [l[, the SM parameters are 



■ 0, 







taken as G F = 1.16638 x lO^GeV - *, a = 1/127.937, m w = 80.385GeV [Uj, m t = 173.2GeV 

0.0405 and |V^Vtd| = 0.0087. The relevant parameters 



12|, m A 



105.6584MeV, \V t * b V ts \ 



of neutral B mesons are collected in table HI Uncertainties of these parameters are not 
considered in the numerical calculation, since the uncertainties induced by these parameters 
on the relevant constraints on the couplings of the MSSM-RPV are far beyond the scope of 
the experimental data. 
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As described in the last section, the contributions of the SM are taken as input, and then 
we calculate the total effects from the sum of the SM and the scalar neutrino contribution 
or the sum of SM and scalar quark contribution of the MSSM-RPV. By demanding the 
total contributions do not exceed the experimental upper and lower bounds, we obtain the 
following constraints on the relevant combinations of the parameters in the MSSM-RPV, 
respectively, 

I y *1^23| < 6 52 x KT 11 , 
-2.29xl0~ 9 < <2.87xl0~ 9 . (17) 

m u i 

The scalar neutrino coupling suffers roughly 2 orders of stronger constraints than the scalar 
quark coupling since there is no helicity suppression in the sneutrino contributions. Due to 
the more stringent experimental limit, our results are several orders of magnitude stronger 



than previous results in the literature 



BE 



6] . If we further assume the mass of the sparticles 
to be several hundreds GeV, a roughly estimate shows that the above combinations of the 
R-parity violating couplings are around 10 -6 and 10~ 4 , respectively, which means that the 
magnitudes of the couplings are not too far away from unity. 

We also give the constraints from the newest experimental upper limits on the B(B® — >■ 
/i + /U~) shown in Eq. (j3J). The corresponding constraints on the sneutrino exchange term 
and the squark exchange term are given below 

|V-AkAa3| < 7 . 8 5 x lo-n, 

i u i 

J-A^3 <117xlQ -8 (lg) 

i m ^ 

IV. SUMMARY 

In conclusion, using the newest experimental data, we have calculated the contributions 
to the B° s{d) -> in the framework of the MSSM-RPV. We gave the constraints on the 

relevant combinations of the parameters in the MSSM-RPV as 

I V- a^3, 6 52 1Q _n 
-2.29 x 10~ 9 < ^^£i|k <2.87xl0" 9 ; (19) 
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E 




< 7.85 x 10 



-11 



E 



< 1.17 x 10" 8 . 



(20) 



Our results are several orders of magnitude stronger than the previous results in the litera- 
ture. 
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